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ABSTRACT 

We present the results of the Hubble Space Telescope i?-Band Parallel Survey (BBPS). It 
covers 0.0370 square degrees and consists of 31 shallow (4-6 orbit), randomly selected high 
latitude HST WFPC2 parallel fields with images taken in both the B (F450W) and I (F814W) 
filters. The goal of this survey is to morphologically classify the galaxies in a homogeneous 
manner and study galaxy properties as a function of type and l?-band magnitude for 18 <; 
bj <^ 23.5 mag. The full sample contains 1800 galaxies, 370 of which are brighter than the 
\ formal statistical completeness limit of bj <^ 23.5 mag. The galaxies are selected from the 

-B-band images and classified using an artificial neural network (ANN) galaxy classifier on the 
higher S/N /-band images. These provide (more) reliable types for / <J 24 mag (or bj <J 26 
mag), since these /-band classifications are less subject to the uncertain redshifted rest-frame 
£f~) ', UV morphology. The ANN classification depends on the shape of the surface brightness profile, 

but not on color. These results are combined with similar (deeper) studies in the Hubble Deep 
Field and the deep WFPC2 field surrounding the radio galaxy 53W002, for which galaxies have 
Q-f been classified to bj <^ 27 mag. 

The galaxy counts for the combined £>-band selected samples show adequate statistics for a 
range 19 <^ bj <^ 27 mag, and are in good agreement with other studies in the flux range where 
they overlap, while showing improved statistics at the bright end. The galaxies are subdivided 
into 3 morphological classes: early- types (E/S0), mid- types (Sabc) and late- types (Sd/Irr), and 
the -B-band counts are presented for each class, as well as the total counts. The faint end of 
the counts is dominated by the irregular galaxies, which have a steep count slope of dlogN / 
dm ~ 0.4. These type dependent counts are compared to models based on local luminosity 
functions which include the effects of the cosmological constant, The whole BBPS sample, 
along with the two deeper fields, is used to delineate the general trends of effective radius and 
(B-L) color as function of both morphological type and apparent magnitude for 18 bj ^ 27 
mag. These properties are discussed in the context of recent redshift surveys. A possible 
explanation for the combined results is given in terms of the effects of £l\ on the evolution of 
the merger rate in a hierarchical scenario. 

Subject headings: Survey — galaxies: fundamental parameters (classification) — galaxies: 
statistics 



1. INTRODUCTION 

Over the last 20 years, the galaxy counts conducted 
in the blue passband showed a remarkable excess of 
faint galaxies relative to model predictions. This ex- 
cess is known as the faint blue galaxy (FBG) prob- 



lem (see reviews by e.g. Koo & Kron 1992; Ellis 
1997). Attempts to model the counts and color- 
distributions of these faint blue galaxies led to the 
conclusion that galaxies were more luminous and 
bluer in the recent past. In order to better under- 
stand the field galaxy population, faint galaxy red- 



1 Based on observations with the NASA/ESA Hubble Space Telescope obtained at the Space Telescope Science Institute, which 
is operated by AURA, Inc., under NASA contract NAS 5-26555. 
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shift surveys were conducted which showed that stan- 
dard luminosity evolution alone cannot account for 
the excess of FBG's (e.g. Broadhurst etal. 1988). 
Many of the faint galaxy spectra show evidence for 
strong star-formation, which led to the conclusion 
that the steep slope of the number counts is produced 
by lower luminosity galaxies undergoing short bursts 
of star- formation (Broadhurst etal. 1988). 

More recently, many ground-based redshift surveys 
have been conducted to faint limits in order to fur- 
ther study this issue. The Canada-France Redshift 
Survey (Lilly et al. 1995a; hereafter CFRS) contains 
591 galaxies brighter than Iab ^ 22.5 mag. They 
have studied the redshifts, emission line strengths 
and ground based photometric properties. With this 
data, Lilly etal. (1995b) showed that the evolution 
in the luminosity function (LF) out to z < 1 was 
greater for the bluer galaxies than for the redder ones 
(presumed to represent late and early morphological 
types, respectively). Recently, this work has been 
extended to include HST based morphology (Brinch- 
mann et al. 1998; Lilly et al. 1998; Schade et al. 1999; 
Le Fevre et al. 1999), as we will use for a larger num- 
ber of field galaxies in this paper. The Canadian 
Network for Observational Cosmology cluster redshift 
survey (Yee, Ellingson k, Carlberg 1996) aimed at 
studying galaxy clusters in the range 0.2 < z < 0.55, 
and was complete for field galaxies with z < 0.3. This 
work has been extended to include more galaxies (Yee 
etal. 2000). For a better measurement of the local 
galaxy luminosity function, the 2dF Galaxy Redshift 
Survey will provide spectra and redshifts for 250,000 
galaxies to a limit of bj = 19.45 mag (Colless etal. 
2001; herafter 2dFGRS). Within several years, the 
Sloan Digital Sky Survey will provide redshifts for 
~ 2 x 10 6 field galaxies to r'=18.0 mag or bj <; 19 
mag (Gunn etal. 1998; Blanton etal. 2001; here- 
after SDSS). These surveys provide very significant 
spectroscopic coverage of relatively nearby galaxies. 
However, as we will show in this paper, the median 
scale- lengths of galaxies at bj ^ 19 mag is r e <j l'.'O, 
and rapidly decreases at fainter magnitudes, so that 
reliable morphological information for faint galaxies 
{bj ^ 19 — 20 mag) over wide fields is beyond the 
capabilities of ground-based facilities and has to be 
done from space. 

The advent of high resolution space-based opti- 
cal imaging opened the door for studying the sub- 
arcsecond properties of many types of astronomi- 
cal objects (e.g. Driver et al. 1995a; Glazebrook 
et al. 1995; Abraham et al. 1996a; Odewahn et al. 
1996; etc). In particular, the Hubble Space Telescope 
(HST) has proven very useful for studying the proper- 
ties of faint, as well as distant galaxies. Most notably 
the Northern Hubble Deep Field (HDF-N, Williams 



et al. 1996) provided deep (I ^ 29 mag) imaging of a 
single field in the UBVI filters at a resolution of bet- 
ter than 071 (actually ~ 0"06 FWHM after drizzling). 
There have been many studies of faint galaxy mor- 
phology from this data set (e. g. Odewahn etal. 1996; 
Abraham etal. 1996b; van den Bergh etal. 1996; 
Driver etal. 1998; Marleau & Simard 1998). The 
HST Medium Deep Survey (MDS) covered a much 
larger area of sky, but to a lesser depth, and took 
advantage of HST's parallel observing mode to im- 
age many fields in mostly V and / (e.g. Ratnatunga, 
Griffiths & Ostrander 1999). This large data set has 
produced many studies (e.g. Driver, Windhorst & 
Griffiths 1995; Abraham et al. 1996a; Im et al. 1996; 
Roche et al. 1996; Im et al. 1999), including one that 
utilized the small fraction of the MDS data that was 
observed in the S-band (Roche etal. 1997). 

The HST data allowed for the morphological prop- 
erties of these faint field galaxies to be studied. One 
of the main discoveries was that most of the FBG ex- 
cess is due to galaxies of irregular or peculiar types 
(Driver et al. 1995a; Driver, Windhorst & Griffiths 
1995; Glazebrook et al. 1995). These authors showed 
that the counts are in excess of the non-evolving LF 
model predictions, regardless of adopted cosmology, 
with the caveat that the local luminosity function has 
to be (arbitrarily) re-normalized by a factor of up to 
2 at a flux level of bj = 18 mag. The justification 
for this re-normalization is typically given as either 
due to: (1) local inhomogeneity (Marzke etal. 1994; 
Zucca etal. 1997); (2) photometric scale errors (Met- 
calfe, Fong & Shanks 1995); (3) rapid recent evolu- 
tion in the galaxy population (Maddox etal. 1990); 
and/or (4) the probability that local surveys are bi- 
ased against late-type low-SB galaxies in selection for 
spectroscopic follow-up. Nevertheless, it remains an 
unsatisfactory situation that an ad hoc correction is 
required to reconcile the local measures of the lu- 
minosity function of galaxies with galaxy counts as 
bright as bj = 19 mags. For more details on this 
issue, see discussions in Driver, Windhorst & Grif- 
fiths (1995) and Marzke etal. (1998). It is therefore 
crucial to have good statistics, with morphology at 
this brighter flux level in the blue filter to attempt 
to explain, or perhaps even rule out, this apparently 
necessary normalization factor. 

Further work on the evolutionary model predictions 
for faint field galaxies has produced conflicting re- 
sults. Without being able to fully cover the literature 
on this topic here, we will try to present the flavor of 
the problem. The bright HDF-N galaxies were used 
to model what the distribution and appearance of 
the fainter HDF-N galaxies should look like assuming 
no-evolution (Bouwens, Broadhurst & Silk 1998a), as 
well as assuming different types/amounts of evolution 
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(Bouwens, Broadhurst & Silk 1998b). Assuming that 
the 30 brightest HDF-N galaxies are a fair represen- 
tation of the thousands of fainter ones, they showed 
that the no-evolution models could not simulate the 
observed properties of the population of the fainter 
galaxies (Bouwens, Broadhurst & Silk 1998a). It was 
then shown that pure luminosity evolution (PLE) 
models fit the galaxy counts, but not the galaxy size 
distribution, and that dwarf augmented populations 
would not fit the counts (see also Babul & Ferguson 
1996; Ferguson & Babul 1998). A model that does fit 
the data was that of "mass-conserving" density evolu- 
tion (Bouwens, Broadhurst & Silk 1998b). These re- 
sults should be viewed with caution, since the HDF-N 
was selected to be devoid of bright objects (Williams 
et al. 1996), which, compounded by the fact that 
galaxies are known to cluster, would logically lead 
to a general deficit of low- to intermediate-redshift 
galaxies in this field. Kauffmann, Chariot & White 
(1996) used the ground-based CFRS colors and red- 
shifts to predict which galaxies are of early-types, and 
then showed that more than the standard amount of 
passive evolution was required by the data. It should 
be noted that the use of colors to determine morphol- 
ogy of distant galaxies was shown to be problematic 
within the same data set, when the HST morpholo- 
gies were added, and no evolution of the space den- 
sity of ellipticals was observed (Schade etal. 1999). 
Another study showed that in an open universe, the 
counts (in the Ubjr^IK bands) and color and redshift 
distributions could be reasonably explained by pure 
luminosity evolution, but that number and luminos- 
ity evolution (NLE) would be required in an = 1 
universe (Pozzetti, Bruzual & Zamorani 1996). In 
contrast, it was shown that PLE models that fit the 
number counts and redshift distributions cannot fit 
the observed (B — K) color distribution regardless 
of chosen cosmology (He & Zhang 1998). This was 
due to the fact that all galaxies with (B — K) > 5.5 
mag must be ellipticals, and the PLE models do not 
match the observed redshift distribution for these red 
galaxies. These authors showed that number lumi- 
nosity evolution better matched all the data that was 
considered (He & Zhang 1998). The luminosity func- 
tions of elliptical galaxies out to z < 1.2, with mor- 
phologies determined from the HST-MDS, were con- 
structed and showed 1 ± 0.5 mag of luminosity evolu- 
tion, and argued against strong number density evo- 
lution (Im etal. 1996). Therefore, it is reasonable 
to say that the issue of the amount or type (luminos- 
ity, density or some combination thereof) of evolution 
is not a solved problem, and is further compounded 
by our poor understanding of the LF normalization 
factor mentioned above. 
Traditionally, galaxies have been classified visually 



by eye. The large number of galaxies that are ob- 
served today necessitated the development of auto- 
mated, computer-based classification methods. Ode- 
wahn etal. (1996) used artificial neural network 
(ANN) classification in a seven-dimensional photo- 
metric parameter space to estimate galaxy stage val- 
ues in the revised Hubble classification system (de 
Vaucouleurs 1959). Abraham etal. (1996a) used lin- 
ear relationships between two parameters, asymme- 
try and central concentration. Marleau and Simard 
(1998) used a 2-dimensional bulge-disk decomposi- 
tion to compare the ratio of bulge-to-total light to 
the types of Abraham etal. (1996a) and of van den 
Bergh etal. (1996). Though the results for a given 
galaxy may differ for any of these methods, the overall 
qualitative results are in reasonably good agreement. 

Where do all these faint galaxy studies leave us? 
Currently, these HST surveys have yielded rather 
poor statistics for galaxy morphology in the mag- 
nitude range bj ~ 18 — 24 mag. Since the cor- 
respnding median scale-lengths of faint galaxies are 
r e ~ l'/O — 0?5 in this magnitude range (see §4.2 and 
figures below), we cannot get reliable galaxy mor- 
phology bj ~ 18 — 24 mag over wide fields from 
the ground. Hence, there remains a need for better 
statistics and good HST morphologies in the mag- 
nitude range of 18 <^ bj 24 mag, where the ear- 
lier surveys have provided little statistics. Since the 
number of galaxies observed at a given magnitude is 
roughly proportional to the area surveyed, and since 
there are more faint galaxies than bright ones, ultra- 
deep surveys like the HDF-N and the field surround- 
ing the weak radio galaxy 53W002 (Odewahn etal. 
1996) have provided important information on galax- 
ies fainter than bj ;> 23 mag. However, the statis- 
tics get increasingly sparse for brighter galaxies. The 
goals, therefore, of this paper are to survey a large 
number of HST fields in the B and I filters over a 
wide area to a lesser depth, and to study the prop- 
erties of these brighter galaxies as a function of ob- 
served .B-band brightness. This will help to fill in this 
new, but essential portion of parameter space and aid 
in the understanding of the faintest galaxies observed 
today. 

The HST observations are described in § 2 and the 
data reduction and catalog generation are outlined 
in § 3. The results of this paper are presented in § 4 
and discussed in § 5 in terms of other recent work in 
the field. We summarize the paper in § 6. For clar- 
ification of nomenclature, we refer to the observed 
radius of a galaxy containing half of the observed 
light as the effective radius (r e ), and this is discussed 
in § 4.2. Since the two deeper fields, HDF-N and 
53W002, are essentially complete past our classifica- 
tion limit at bj £ 27 mag (Odewahn etal. 1996), we 
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shall group them together and refer to them as the 
.B-band Deep Survey (BBDS), while the 31 shallower 
fields will be referred to as the -B-band Parallel Sur- 
vey (BBPS). We adopt H Q = 65 km s^ 1 Mpc~ l and 
a flat [VI m = 0.3, VIa = 0.7) cosmology throughout, 
except where otherwise noted. 

2. OBSERVATIONS 

All data for this survey were taken with the Hubble 
Space Telescope using the Wide-Field and Planetary 
Camera 2 (WFPC2) in parallel mode during Cycles 
6-7. Fields were randomly selected with the criteria 
that they be at high Galactic latitude (| b 11 | > 30°) 
and that they contain no bright SAO stars, RC3 
galaxies or known Abell galaxy clusters. The data set 
consists of 31 fields, all selected to have low Galactic 
extinction, Ab < 0.4 mag (Burstein & Heiles 1982). 
The data are summarized in Table 1. Note that fields 
bb025 and bb026 were not used in the final analysis 
because they contained single /-band exposures that 
were inadequate for this survey. 

For each field, images were taken in both F450W 
(the WFPC2 5-band) and F814W (the WFPC2 I- 
band) with a longer exposure time in B due to the 
approximate 20% lower sensitivity at that wavelength 
(see Table 8 of Holtzman et al. 1995). Combined with 
the expected relative colors of the faint field galaxies 
compared to the Zodiacal sky, the total survey ef- 
ficiency in F814W is ^ 2x better than in F450W. 
Therefore, we took 2-4 orbits in B and two orbits in 
/. The reason for using both filters is to select the 
galaxies from the S-band images, where the HST res- 
olution improves deblending over ground-based data, 
and therefore provides the most reliable object se- 
lection and photometry. The I-band images (which 
corresponds to the rest-frame -B-band) have signifi- 
cantly higher S/N, and therefore provide the best pos- 
sible source of morphological classifications, as well as 
provide color information over a wider color baseline 
than the standard V-I color (F606W-F814W) that 
has traditionally been used in most WFPC2 studies. 
The two sets of images coupled with a future redshift 
survey will also allow us to do a study of the rest- 
frame color gradients in the brighter galaxies. The 
total area of the BBPS survey, including the PC ex- 
posures, is approximately 0.0370 square degrees. The 
WF(PC) data have a / /0996(0 / /0455) spatial sam- 
pling and a limiting magnitude of bj <^ 23.5 mag (the 
approximate 90% completeness limit for compact ob- 
jects), as discussed in § 4.1. 

3. DATA REDUCTION 

3.1. Image Stacking 

All images were registered using the centroids of 
manually selected compact and bright objects, with 



the centroids being determined from the IRAF task 
imcentroid. Since the detectors are fixed relative to 
each other, shifts were determined for each of the 4 
WFPC2 CCDs, and then averaged together to get 
the optimal offset for each set of exposures in a given 
field. These offsets were rounded off to whole num- 
bers and only integer shifts were applied, so as not 
to introduce additional numerical noise into the data 
(c/. Ratnatunga, Griffiths Sz Ostrander 1999). Due 
to the nature of parallel observing with HST, a set 
dither pattern cannot be chosen by the parallel ob- 
servers, and so the method employed here produces 
the most self-consistent data set. The appropriate 
shifts were then applied taking into account the dif- 
ferent chip orientations and the finer sampling of the 
PC. 

Since most fields were two orbit cases, it was nec- 
essary to improve upon existing cosmic ray (CR) re- 
jection algorithms, which were optimized for use with 
N > 5 images (e.g. Windhorst, Franklin & Neuschae- 
fer 1994). All image stacks were created using a 
customized IDL routine which was specially devel- 
oped for this project. This routine was also used for 
the low light-level images in the F410M filter of Pas- 
carelle, Windhorst & Keel (1998). Images through 
different filters were handled separately to ensure that 
our final morphological classifications were as color- 
independent as possible. The reason for developing 
this new routine was the need to accurately reject 
CRs over a few (2-4) independent exposures, while 
assuring at high reliability that the science image it- 
self would not be corrupted by the algorithm. 

For N registered images, the IDL routine performs 
the following operations at each pixel (x,y) location 
separately. It creates a list of N pixel values which is 
then sorted from the lowest to the highest value. The 
following Poisson noise model based on the known 
CCD characteristics is then used to determine which 
pixel values should and should not be included in the 
average: 

<Tx, y = \]DN x , y * g + RN 2 + DK * t/g (1) 

where DN XiV is the number of ADU in pixel (x,y), 
g = 7.0 e~/ADU is the WFPC2 detector gain, 
RN = 5.3e~ is the read-noise, DK = 0.0033e~ /sec 
is the dark current rate and t is the exposure time 
in seconds. Starting with the lowest pixel value, each 
successive value is checked to see if it is within 2.5 a 
of the current average value. This way, higher pixel 
values that are likely due to cosmic rays are rejected 
(at the 2.5 a level). This process is then repeated 
for each pixel. This rejection algorithm will fail for 
the two orbit cases when a pixel is hit by a CR in 
both images, which we measured to occur about 0.3% 
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of the time, or for about ~ 2000 pixels that were 
affected by CR's in both full-orbit WFPC2 CCD ex- 
posures. This number is in excellent agreement with 
the values given in the WFPC2 Instrument Hand- 
book 2 (Biretta et al. 2000) when extrapolated to our 
longer exposure times. These left-over CR's cause 
problems with automated image extraction, because 
they would be counted as objects and possibly clas- 
sified as faint galaxies, which would contaminate the 
sample. 

In order to deal with this problem, we checked for 
cases of "double hits" by subtracting a 3x3 pixel 
median-smoothed image from each original image, 
and compared this difference image to the difference 
image for the other orbital exposures. This way, all 
"double hits" were located and interpolated over by 
substituting the median of the 8 surrounding pixels. 
The result was a much cleaner looking image, but this 
process also unavoidably generated holes in the cen- 
ters of bright, centrally-concentrated objects (stars 
and early- type galaxies). Since these holes would 
obviously contaminate the photometry of the bright 
galaxies we are interested in, we opted to use the 
cleaner but non-photometric images for the object se- 
lection only. Once the object positions were found, all 
analysis was performed on the images that still con- 
tained the "double hit" CR's, but that had objects 
that were uncorruptedin their central regions. 

3.2. Object Extraction 

All image detection was done using the SExtractor 
version 1.0a (Bertin & Arnouts 1996) object finding 
software package. Each WFPC2 WFC (or PC) image 
is convolved with a 7x7 (or 9x9) convolution mask 
of a Gaussian kernel with a FWHM of 3 (or 5) pixels. 
Then all objects with at least 8 contiguous pixels that 
are 2.5 a above the local sky-background (as deter- 
mined by SExtractor) are extracted, and their loca- 
tion, size, and magnitude are written to a file. The 
mean isophotal 2.5a detection limits for each field are 
listed in Table 1. Averaged over all fields, this sur- 
face brightness ("SB") detection limit corresponds to 
SB(6j)~ 24.2 ±0.2 (22.8 ±0.2) mag arcsec" 2 for the 
WF (PC) detectors. The la SB sensitivity limits for 
53W002 field of Windhorst, Pascarelle & Keel (1998) 
are SB(6j)~27.5 mag arcsec -2 and SB(I)~26.7 mag 
arcsec -2 , and these numbers are about ~1.5 mag 
fainter for the HDF-N (Odewahn etal. 1996). 

For each detected object, a smaller sub-image is 
then cut whose size is determined by the SExtrac- 
tor image parameters. The subsequent image anal- 

2 Updates to the WFPC2 Instrument Handbook can be found at 

3 Kron (1980) gives the transformation from standard to photo, 
Holtzman et al. (1985) Table 10, we derive B - F450W = 0.23(B 
the same as the photographic bj system. 



ysis is done on these individual object images (i.e. 
all surface photometry and morphological classifica- 
tion) , using the MORPHO package of Odewahn et al. 
(1996, 1997). 

Once the sub-images are extracted, all objects are 
visually inspected to determine whether they are true 
detections, as well as to determine if SExtractor over- 
did the deblending of neighboring objects. The lat- 
ter was especially a problem for the brighter spiral 
and irregular galaxies that had several peaks in their 
2-dimensional brightness distributions, and is a gen- 
eral problem for other object finding algorithms as 
well (see e.g. Valdes 1982; Neuschaefer &: Windhorst 
1995). The number of objects that had to be re- 
extracted by hand was relatively small, about 10 per 
field. The number of separate objects that were close 
enough to overlap was also small, but SExtractor gen- 
erally did an excellent job of deblending these, al- 
though we emphasize that any deblending process is, 
by definition, somewhat arbitrary. The photometry 
of these blended objects is addressed below. 

3.3. Surface Photometry 

Photometric zero-points were taken from Tables 9 
and 10 of Holtzman etal. (1995). In the visual- 
red, these have not changed much during the life- 
time of WFPC2 (Biretta etal. 2000). We used the 
"synthetic" WFPC2 zero-points (Vega-based), which 
means that the F450W magnitudes are equivalent to 
bj 3 , and the F814W magnitudes are the same as the 
ground-based / magnitudes. This is to allow optimal 
comparison to previous ground-based galaxy l?-band 
counts, which were primarily done in the bj filter (or 
a filter that easily transforms to bj). For complete- 
ness, we used zeropoints of ZPp^ow = 21.929 mag 
and ZPpsuw = 21.676 mag for bj and I, respectively, 
for a count of 1 DN per second of exposure time. 

In order to accurately derive light-profiles and com- 
pute total magnitudes, a good determination of the 
sky-background must be made. For each field and 
filter, sky-values are determined for each of the four 
WFPC2 CCD's using the following technique. The 
CCD is divided into 16 squares of 101x101 pixels 
each. Next, a sky- value and sky-sigma are deter- 
mined for each square by fitting a "tilted plane" (c/. 
Neuschaefer & Windhorst 1995) to that area, after 
rejecting all real objects that were found at the > 2a 
above the local sky level. The reason for not fitting a 
more complicated surface is that the images are suf- 
ficiently flat-fielded over ~ 100 pixel scales, so that 
mostly linear sky gradients are left, if any. Resid- 
http: / / www.stsci.edu / instruments / wfpc2 /wfpc2_top.html 
graphic B magnitudes as B — bj = 0.23(B — V) and, using 
V), which shows that the synthetic WFPC2 F450W system is 



6 



£>-Band Parallel Survey 



ual sky gradients could be due to, e.g., the telescope 
pointing too close to the Earth's limb, which is not 
always controllable when scheduling parallel observa- 
tions. This could produce a simple mono-directional 
linear gradient in the WFPC2 sky, although this was 
not often observed. The median of the 16 sky-box 
values was then used to determine the final sky- value 
and sky-sigma. These values were found to agree with 
the global BACKGROUND and RMS values deter- 
mined by SExtractor to within 2%. 

Since the Zodiacal sky at moderate to high Eclip- 
tic latitudes is about SB(6j)~23. 8-24.1 mag arcsec -2 
and SB(I)«22. 15-22.45 mag arcsec -2 (Windhorst, 
Keel &: Pascarelle 1998), a 2% error in the sky- 
subtraction becomes the dominant factor at surface 
brightness levels SB(6j)~28.0-28.3 mag arcsec -2 or 
SB(I)«26.4-26.7 mag arcsec -2 , at which level this 
sky-subtraction error equals ~100% of the galaxy SB- 
profile signal. In the noisier BBPS images, we cannot 
push the galaxy SB-profiles this faint, but for the two 
BBDS fields this is one reason why we do not push 
the galaxy classifications (see §4.1) fainter than total 
fluxes of bj = 27.0 mag (or I = 25.5 mag). This is 
because at these total flux levels, the median galaxy 
scale-lengths are r e = C/.'2 — C/.'3 (see §4.2), so that the 
sky-subtraction errors in the average SB-value out to 
r e are of order 10-20%. This is the maximum accept- 
able error in SB-profile parameters to allow reliable 
ANN classifications, which are based on the galaxy 
SB-profile and other parameters (see §3.5). Fainter 
than bj ;> 27 mag, ANN galaxy classifications are 
mainly unreliable due to a lack of a good training 
set (Odewahn et al. 1996), which would be neces- 
sarily based on the abilities to assign eyeball classi- 
fications as well as accurately measure the light pro- 
files of these galaxies, both of which become increas- 
ingly suspect fainter than this flux level, for reasons 
stated above. The issue of uncertain resdshifted UV- 
morphology could also play a role here and is further 
discussed in § 3.7.2. Given that bj ~ 27.0 mag is 
also the limit due to the sky background subtraction 
(at these wavelengths), there is little hope of reliable 
classification for fainter galaxies with HST, even with 
improved resolution and/or deeper images. 

All surface photometry was performed on each sub- 
image, using these sky-values, with the MORPHO 
package as described in Odewahn etal. (1996, 1997). 
Since each sub-image to be analyzed has to contain a 
single object, the analysis on crowded areas was done 
as follows. If two objects were too close together, two 
sub- images were generated, each with one of the two 
objects replaced by the local sky value, as interpo- 
lated from surrounding unaffected sky areas. The size 
and shape of the "blanked" out region was an ellipse 
based on the SExtractor values for position, size, axis 



ratio and position angle. This way, when the isopho- 
tal ellipses were fit to determine the light-profile, the 
target object was minimally affected by its neighbor. 
Since the total magnitude is determined from the in- 
tegrated light-profile, this was the most reliable way 
to handle proximate objects, while reducing the ef- 
fects of crowding on the other measured object pa- 
rameters. Note that because the ANN classification 
is based on individual object light-profile based pa- 
rameters (§ 3.5), and simultaneous light-profile fit- 
ting on more than one object is very difficult to do 
(Schmidtke etal. 1997), this seemed to be the most 
robust way to get an accurate ANN class for every 
object, even if close neighbors were present. 

3.4. Catalog Generation 

Once the images in both the B and I filters are 
processed and cataloged, the B and I object lists are 
cross-matched by pixel location. For each object in 
the .B-band catalog, the I-band catalog is searched 
using a search radius of 3 effective radii. The ef- 
fective radius r e used here is the half-light radius as 
measured from the -B-band image. With this search 
criteria, about 70% of the objects in the .B-band cat- 
alog have /-band counterparts, while more than 95% 
of the I-band catalog have -B-band counterparts. A 
manual check of the images showed that almost all 
the unmatched .B-band catalog objects were not real 
galaxies, but likely image defects, noise spikes or re- 
maining cosmic rays left in the .B-band images. We 
thresholded at 2.5a so as not to exclude any real ob- 
jects, at the expense of making false detections which 
would get thrown out in cross-matching. Since the 
I-band images are of higher quality and go deeper, 
despite their generally shorter exposures, it is doubt- 
ful that very blue objects were excluded through this 
selection process. However, this matching procedure 
could select against extremely red objects, which we 
must bear in mind during the interpretation (§4 and 
5). 

3.5. Morphological Galaxy Classification 

Galaxy types are assigned to all objects using an 
automated artificial neural network (ANN) pattern 
classifier that is based on the shape of the object light- 
profile. Galaxy types, or hereafter T-types, were des- 
ignated on the 16 step revised Hubble system (de 
Vaucouleurs 1959). This classification scheme was 
chosen because T-types have been shown to corre- 
late well with physical properties in nearby galaxies 
(Roberts & Haynes 1994) and the available data used 
to produce our galaxy count models are segregated 
in the same way. Since the typical machine classifica- 
tion error in type is 2-3 steps (Odewahn etal. 1996), 
types were rebinned into three broad intervals, that 
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each span a much larger range than the ANN classifi- 
cation error (E/SO, Sabc, Sd/Irr). This classification 
error results from noise in the profile-parameters of 
the faint objects (as discussed in § 3.3), and from 
inconsistencies between experienced classifiers, who 
visually classified a limited number of galaxies in the 
"training set" that is used to train the ANN be- 
fore it is run on a much larger sample. To the ex- 
tent possible, the personal offsets and biases between 
these expert classifiers were removed (cf. Odewahn 
etal. 1996; Odewahn etal. 2002) before constructing 
the final training set, thereby reducing systematics 
in the ANN classifications (although not yet proving 
that these classifications are necessarily correct). The 
ANN's used here are the rest-frame networks devel- 
oped in Odewahn etal. (1996), which were trained 
on the HDF-N galaxies that had published redshifts 
at that time. The redshifts are required to know the 
rest-frame wavelength in which each galaxy in the 
UBVI training set was actually observed. T-types 
are in the range — 7 < T < —0.5 for objects desig- 
nated as "E/SO," -0.5 < T < 5.5 as "Sabc", and 
5.5 < T < 12 as "Sd/Irr". 

We note that the ANN's of Odewahn etal. (1996) 
always produces an answer, i.e., based on the actual 
galaxy SB-profile parameters, it forces a classification 
in one of these three T-type intervals, no matter how 
faint or how low-SB the object is. That is, we delib- 
erately sacrifice classification accuracy to get 100% 
completeness in the classifications, so that no objects 
remain unclassified. The corollary of this is, of course, 
that beyond a certain flux and SB limit, the classifi- 
cations become increasingly unreliable due to the lack 
of an adequate training set. The lack of an adequate 
training set at these faint flux levels is due to a combi- 
nation of insufficient S /N and resolution for a human 
classifier to assign a reliable type, as well as a poor 
understanding of the rest-frame UV morphology for 
the higher redshift objects where this is a potential 
issue. As discussed in §3.3, we believe that this limit 
occurs in the two BBDS fields at total fluxes of ap- 
proximately bj = 27.0 mag (I = 25.5 mag). We will 
show in §4.1 that the shallower BBPS fields produce 
consistent classifications close to their completeness 
limit of bj ~ 23.5 mag in the flux range of overlap 
with the two BBDS fields (22.5 £ bj £ 24.0 mag), 
where the classification of the latter objects are ro- 
bust, since they have much higher S/N at the same 
flux level. 

Using the suitably trained neural networks, galaxy 
types were generated from both the -B-band and I— 
band images, as described below. However, for a va- 
riety of reasons, all types used in this paper are based 
on those ANN-values assigned from the I-band im- 
ages. First, the I-band images are superior to the B- 



band ones due to the larger overall sensitivity and the 
much darker Zodiacal sky background at that wave- 
length, even for the bluest detectable field galaxies. 
The second reason is bandpass shifting. The median 
redshift of the fainter galaxies at bj < 25 mag is 
z ~ 0.6 - 0.7 (Koo & Kron 1992; Ellis 1997). There- 
fore, the light observed in the I-band was emitted in 
approximately the -B-band where morphology of lo- 
cal galaxies has been well studied, and a good train- 
ing set for the ANN exists. Analogously, the light 
observed in the -B-band is emitted in the rest frame 
near-UV, where systematic studies of nearby objects 
have only recently been done ( cf. Giavalisco et al. 
1996; Burg etal. 1997; Kuchinski et al. 2001; Mar- 
cumetal. 2001; Windhorst et al. 2002). The issue 
of uncertain rest-frame UV morphology and its effect 
on the galaxy counts as a function of type is further 
discussed in § 3.7.2. For the BBPS, the use of I-band 
classifications allows us to take advantage of bandpass 
shifting instead of being hampered by it. Thirdly, the 
effects of the k-corrections are smaller at this longer 
wavelength. Finally, the HST images are closer to 
being properly sampled in the I-band, compared to 
the -B-band, allowing for better classifications. 

We classified all these galaxies from their I-band 
images, but using an ANN defined at a shorter wave- 
length which is closest to the rest-frame wavelength 
sampled by the I-band for that particular object. 
These rest-frame ANN's were constructed from the 
HST UBVI images in the HDF-N for objects with 
spectroscopic redshifts (z <; 1). The details of this 
classification method and the use of rest-frame ANN's 
are given in Odewahn et al. (1996, 1997) and refer- 
ences therein. For this project, seven parameters are 
used as input to the ANN which then produces a 
type. These parameters are all based on the object's 
azimuthally averaged radial surface- brightness pro- 
file and are: 

1. the SB at the radius containing 25% of the total 
light; 

2. the SB at the radius containing 75% of the total 
light; 

3. the mean SB within the effective (half-light) ra- 
dius; 

4. the mean SB within an isophotal radius 
(SB (I) = 24.5mag arcsec" 2 ); 

5. the slope of a linear fit to the profile in r 1 / 4 
space; 

6. the intercept of a linear fit to the profile in r 1 / 4 
space; and 

7. the axial ratio (b/a) of the outer isophote. 
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It is important to note that in the ANN classifica- 
tions, color is not used at all and that the effective 
radius is not used directly. Fig. 1 shows that this 
indirect use of effective radius does result in a loose 
correlation of r e with type in any given I-band mag- 
nitude slice. It is also important to reiterate that no 
bulge-to-disk decomposition or model profile-fitting 
is used. This particular choice of the 7 ANN param- 
eters is somewhat arbitrary. Different sets of struc- 
tural parameters were tested (Odewahn et al. 1996) 
and produced similar classifications to within the er- 
rors quoted above. 

As discussed in Odewahn etal. (1996), the classifi- 
cation limit in the two BBDS fields is bj ^ 27 mag (or 
/ ^ 25.5 mag). At fainter flux levels, human classifi- 
cation becomes increasingly unreliable. As illustrated 
in § 4 (Figures 3-7), it is interesting to notice that the 
ANN classifications appear reasonable, even a mag- 
nitude or so fainter than the visual classification limit 
of the training sets. 

These classifications are reasonable because the pa- 
rameters that are not used in the classifications (color 
and effective radius) follow the type-dependent trends 
for galaxies at cosmological distances that one would 
expect from the known properties of local galaxies 
(see discussions in § 4.2 and § 4.3). However, this is 
not proof that these classifications are therefore cor- 
rect, and two things are required to verify these clas- 
sifications beyond the currently stated limits: higher 
resolution images of these galaxies and a better mid- 
UV training set for the ANN (see § 3.7.2). For the 
brighter BBPS sample, the /?-band data is too noisy 
near its completeness limit to believe the classifica- 
tions from the S-band images alone. This is another 
reason for classifying in the /-band images, which go 
£3 0.5 — 1.0 magnitudes deeper even for blue galax- 
ies, and therefore have a higher signal-to-noise ratio. 
Therefore, all BBPS fields have reliable classifications 
to their respective S-band detection limits listed in 
Table 1, and the two BBDS fields have reliable clas- 
sifications to bj <^ 27 mag (ie. / <; 25 mag). 

3.6. Star- Galaxy Separation 

The issue of star-galaxy separation is an important 
one, especially in a survey of this nature where the 
goal is to count different types of objects, including 
the most compact ones. The star- galaxy separation 
for the BBPS was initially performed using a linear 
division in a single 2-D photometric parameter space. 
The effective radius, r e forms the Y axis. The X 
axis is formed with the (MORPHO) "C42" param- 
eter. This is a concentration index formed by the 
ratio of the size of the "total aperture" to the 50% 
flux aperture (i.e. the number of effective radii inside 

4 Fortran code to compute Bahcall & Soneira Galaxy Model is 



the total aperture). Both apertures have an elliptical 
shape defined by the median ellipticity of the ellipse 
fits to the lower SB isophotes. The determination of 
the "total aperture" is a complex procedure, which 
finds the optimally sized aperture based on the local 
signal-to-noise ratio of the sky. An object is classi- 
fied as a star if C42 > 60.0, or if it lies below the 
line r e = 0.010 * C42 - 0.095. These values were 
determined using the F814W data from the HDF- 
N, the HDF-N flanking fields, the 53W002 field and 
the 6 earlier MDS fields from Driver, Windhorst & 
Griffiths (1995). They were plotted for several mag- 
nitude cuts using symbols to indicate the morpholog- 
ical types assigned in our past visual classification 
work. Bright objects with clear diffraction spikes 
were always visually classed as stars. These could 
erroneously include a few QSO's, if these were not 
otherwise recognized, i.e., from their available spec- 
tra (e.g., Pascarelle etal. 1996) or weak radio fluxes 
(Windhorst, Keel & Pascarelle 1998). By inspect- 
ing many such plots, the best line that divided the 
stars and the E/S0 types was determined. At bright 
magnitudes (usually I < 22 mag for most BBPS 
fields), the parameter space always showed a fairly 
substantial segregation between stars and even com- 
pact, high SB galaxies (likely E/S0). Of course, as we 
go fainter, this division blurs due to increasing noise 
in the extracted r e and C42 measurements. This sep- 
aration limit is similar to the one found by Mendez & 
Guzman (1998), who used a different method. Also, 
a visual inspection of the 40 objects with bj < 21 mag 
and not classified as stars showed that this method 
failed for 12 bright (bj < 19 mag) and saturated stars. 
These 12 stars were removed from the the final galaxy 
counts. 

To demonstrate the success of this star-galaxy sep- 
aration method, Fig. 1 shows the /-band magnitude 
versus effective radius relation for the BBPS objects, 
including those objects that were classified as stars. 
It is clear that to the BBPS /-band compact-object 
detection limit (/ ~ 23.75 mag), the stars easily sep- 
arate out from most of the galaxies, except for an 
interesting group of objects, most of which were ini- 
tially classified as "elliptical" galaxies, that are small 
(log(r e ") ~ — 1) and near the faint end of the sample 
(/ ~ 22.0 mag - these objects are plotted as pluses in 
Fig. 1). This raises the obvious and important ques- 
tion of whether or not these objects are misclassified 
stars. 

In order to test this issue, the Galaxy Model 4 of 
Bahcall & Soneira (1981) was used to predict the 
star-counts in these fields. The model was run for 
the (I 11 ,b n ) coordinates of each of the 29 telescope 
pointings and the resulting predicted star counts were 

available at http://www.ias.edu/~jnb 
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averaged together. Fig. 2 shows the differential num- 
ber counts in the i?-band for all non-stellar objects, 
for the objects classified as stars, for the E/SO galax- 
ies only, as well as for the average Bahcall and Soneira 
model. The first issue to notice is that the average 
star count model matches the observed star counts 
very well. When compared to the star count model, 
it is clear that our first pass of the star counts (plotted 
as upwards arrows in Fig. 1) turn over at bj > 23.5 
mag due to a significant fraction of misclassifications 
and incompleteness, even before the galaxy counts 
noticeably turnover at bj > 24.5 mag. It is there- 
fore highly probable that some fainter stars have 
been misclassified as galaxies. These would have 
most likely been misclassified as centrally concen- 
trated (and therefore elliptical) galaxies, as seen in 
Fig. 1. There is also a color selection effect, in that 
the fainter stars are expected to have redder colors 
and therefore may have escaped the blue selection 
limit. Note that, in order for an object to appear in 
our catalogs, and hence in Figs. 1-7, it had to be seen 
in both the B and I filters. The majority of the stars 
from the two deeper fields with I > 22 mag are too 
red [(B — I ) > 4 mag ] to have been selected in the 
BBPS sample. This is less of a problem for the galax- 
ies because, as Fig. 7 shows, at I > 22 mag galaxies 
have (B—I) colors in the range of 1.0-2.5 mag. For- 
tunately, at this brightness level and at fainter flux 
levels, the galaxy counts are much higher than the 
predicted star counts (see Fig. 2), even for the early- 
type galaxies which have the shallowest counts slope 
(see also Fig. 3). Therefore, the likely contamina- 
tion of the faint galaxy counts by misclassified stars is 
minimal and not significant. The few missing bright 
(bj <^ 18 mag) stars in Fig. 2 are due to the satura- 
tion problem mentioned above. These missing bright 
stars were obviously detected, but due to saturation, 
no reliable fluxes could be measured, so they are not 
plotted in Fig. 2. 

Since the star-counts should go as deep as the to- 
tal object counts, and Fig. 2 shows that our first-pass 
star counts (shown as upwards arrows) turns over be- 
fore the galaxy counts, there must be some number of 
fainter stars that were not classified as such. There- 
fore, a second step in star-galaxy separation was ap- 
plied to rectify this situation. The model star counts 
were used as a reference point to update the stellar 
classification for bj > 23 mag. Therefore, the result- 
ing star counts can not and must not be used fainter 
than bj > 23 mag to test any star count models. This 
is not a problem here, since the goal of this study is 
to properly count the galaxies in order to compare 
to galaxy count models as a function of galaxy type. 
A diagonal line was drawn in the magnitude-log(r e ) 
parameter space, where both quantities are measured 



in the /-band. Note that Fig. 1 shows the effective 
radius (r e ) as an average over the values measured 
in the B and I bands, while all star-galaxy separa- 
tion (and galaxy classification) was done exclusively 
on the I-band images. All objects to the left of this 
line were declared stars and all objects to the right 
were declared galaxies. This step did not change the 
classification of any previously determined stars. The 
process was iterated by trial and error until the ob- 
served star counts matched the predicted star counts 
up to the completeness limit of the total counts. The 
final star-galaxy separation line that was used is given 
by I = -20.313 x log(r e ) + 6.188. The objects that 
were re-classified as stars are shown as superposed 
pluses in Fig. 1, and it can be seen that the majority 
of them were originally classified as ellipticals. 

However, for any reasonable application of this 
second pass at star-galaxy separation, the new star 
counts were always several sigma above the Bahcall 
and Soneira models (at the faint end) and so war- 
ranted further investigation. The colors of these faint 
centrally concentrated objects were examined and 
a significant number of them were grouped around 
{(B — I)) = 1.44 ± 0.42 mag and the majority of 
these came from 2 HST parallel fields in the Virgo 
Cluster, bb016 and bb018. Using the NASA/IPAC 
Extragalactic Database (NED), it was realized that 
field bb016 is only 8'.2 away from M87. This well 
studied elliptical galaxy is known to have a large and 
extended globular cluster population (Harris 1986). 
If we designate all the compact objects in this field 
with the colors mentioned above as globular clusters 
in Virgo (loosely associated with M87), and compare 
the number that we find (N=19) to the surface den- 
sity of these objects in Harris (1986), we get very 
reasonable agreement. The apparent magnitude of 
these objects, combined with the M87 distance, also 
argues that these have the expected luminosities of 
Virgo globular clusters. Inspection of the location of 
the bb018 field shows that it is situated between two 
bright Virgo ellipticals M84 (8'.2 away) and M86 (9f9 
away). It is therefore highly likely that the bb018 
star-like objects (N=18) in this color range are glob- 
ular clusters in the (overlapping) halos of both of 
these galaxies. The final corrected star counts, shown 
in Fig. 2, have these globular cluster candidates re- 
moved, and show that the updated star-galaxy sep- 
aration has no noticeable effect on the total galaxy 
counts, but do impact on the counts of the faintest 
ellipticals, albeit minimally. In fact, given what is 
known about faint star counts (Flynn, Gould & Bah- 
call 1996), even if all faint stars were mis-classified 
as early- type galaxies, the E/SO counts would be at 
most 10% too high in this flux range. 

We emphasize that star-galaxy separation with 
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HST is non-trivial, even at the relatively bright limit 
of I 22 mag, because faint galaxies generally have 
small angular sizes (see §4.2 and Fig. 6 below), and 
because objects other than stars and galaxies can 
sneak into the HST samples. 

3.7. Galaxy Counts as a Function of Type 

Objects classified as galaxies were counted in half- 
magnitude wide bins. The counts for each field were 
combined in the following way. The completeness 
limits in each field is different for both observational 
(mainly varying exposure times) and statistical rea- 
sons (i.e., cosmic variance). First, the completeness 
limit for each field is determined from the flux level 
at which the total counts in that field turn over (cf. 
Neuschaefer & Windhorst 1995). All counts beyond 
the completeness limit are given a weight of zero while 
all counts brighter than this limit are given a weight of 
one. Therefore, no field contributes to any magnitude 
bin (in the combined counts) in which it is less than 
90% complete. The center of the faintest bin that 
was used from each field is listed in the last column 
of Table 1. This ensures that the fainter points are 
properly handled given the different exposure times. 
These weights are then used to keep track of the area 
that contributes to the counts in each bin. This re- 
sults in some of the fainter bins having a smaller ef- 
fective area than the brighter bins, as indicated in 
Table 2. This is not a problem because there are 
many more faint objects and hence the statistics in 
the faint bins are not compromised. Below we discuss 
two possible sources of contamination, other than the 
stars already mentioned in § 3.6. 

3.7.1. Local Large Scale Structure 

Some of the parallel fields had coordinates that 
placed them in or near the edge of the Coma or 
Virgo superclusters (as indicated in Table 1). In or- 
der to test whether there was any excess of bright 
galaxies from these superclusters, the counts were 
compared for the cluster and non-cluster fields. For 
fainter flux levels, the counts from the cluster fields 
were within the 3a field-to-field variance of the mean 
of the non-cluster field counts. The only statisti- 
cally significant contaminants were two galaxies from 
bb019 (bj = 17.08, 17.22 mag) and one from bbOOl 
(bj = 17.97 mag), which are in the Coma superclus- 
ter. This was the primary HST target which may, 
therefore, have biased the bright galaxy sample in 
the WFPC2 parallels we are studying here. This a 
valid and well known concern in HST parallels (cf. 
Casertano etal. 1995), and the obvious remedy is to 
excise these excess bright galaxies from the sample. 
These three galaxies were therefore omitted from the 
final counts because they are not representative of 



field galaxies in general, given that the primary HST 
target was biased towards these clusters. 

Also the brightest galaxy in the survey (bj = 16.99 
mag) was excluded from the counts because it is 
part of the NGC 383 group, which was also the pri- 
mary target of that particular HST acquisition. This 
galaxy is an elliptical that appears to be interacting 
with a fainter spiral galaxy, and it was only noticed 
because the ANN assigned a much later type. A care- 
ful inspection shows that this galaxy contains a signif- 
icant disk structure with low level spiral arms which 
explains why it was incorrectly classified. Therefore, 
there was only one true "field" galaxy found in this 
survey brighter than bj <J 19 mag. 

To let the readers judge our line of reasoning for 
themselves, we plot the three bright data points based 
on these four galaxies in the total counts (upper left 
panel in Fig. 3), but plot their detected surface den- 
sity as upper limits, since the true field galaxy count 
at these flux levels (16.5 <J bj <^ 18.0 mag) must be 
lower than the count that we observe in these few 
not-so-random parallel fields. This is important for 
the discussion of the counts in § 4 and 5. 

3.7.2. Misclassification Trends due to the Uncertain 
Rest-frame UV 

The mid-UV (2500-2900A) images of nearby galax- 
ies available thus far (Windhorst et al. 2002) suggest 
that it is more likely to misclassify true ellipticals or 
early-type spirals in the rest-frame mid-UV as late 
type galaxies than the other way around. This is be- 
cause truly late-type galaxies are dominated by young 
and hot stars in filters from the mid-UV to the red, 
and so have to first order the same morphology and 
a rather small morphological k-correction. However, 
early-type galaxies (ellipticals and early-type spirals) 
can, although do not have to look dramatically dif- 
ferent when one goes from the rest-frame mid-UV to 
the optical-red part of the spectrum. Hence, misclas- 
sifications due to the morphological k-correction will 
more likely move more ellipticals and early-type spi- 
rals into the late type/irregular category, rather than 
vice versa. Since all galaxies in the present study 
were classified from /-band images, this would only 
be relevant for galaxies with redshifts of z ;> 1.8, 
which corresponds to the fainter part of our sample, 
as can be seen in the morphologically segregated red- 
shift distributions of Driver etal. (1998). These red- 
shift distributions show that significant numbers of 
z ^ 1.8 galaxies are only present in fainter sample 
where the classifications become increasingly unreli- 
able (bj ^ 27.0 mag) for reasons discussed in § 3.5. 
While this can perhaps explain part of the excess late- 
type/irregular galaxies at faint magnitudes as dis- 
cussed below, it cannot explain all of the FBG excess, 
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and it certainly cannot explain the apparent excess of 
early-mid types at bj ^ 24 mag relative to the best 
available models, as discussed below. 

4. RESULTS 

Here we present our results for the ~ 1800 BBPS 
galaxies function of measured i?-band bright- 

ness and ANN morphological type. The statistically 
complete sample is composed of ~ 370 galaxies with 
19 £ bj & 23.5 mag. 

4.1. The Morphological Galaxy Counts 

The fundamental question we will address in this 
section is that of the galaxy counts as a function of 
morphological type at fluxes between the ultra-deep 
HST fields and what can be done from the ground. 
We will show that, in order to use galaxy counts as 
a function of morphological type to address issues of 
galaxy evolution, this intermediate flux range is cru- 
cial. 

In Fig. 3, we show the BBPS galaxy counts as a 
function of morphological type. Also included for 
comparison are the deeper counts from the two BBDS 
fields (Odewahn etal. 1996). At the bright end 
(bj <J 18 mag), we include as upper limits the four 
galaxies, discussed in § 3.7.1, that were excluded from 
the sample due to their high probability of not being 
representative of field galaxies. Also plotted are the 
combined ground-based counts from the Millenium 
Galaxy Catalog (MGC; Liske etal. 2002), which is 
a wide-field l?-band CCD survey covering 30 square 
degrees, which we have analyzed using the same soft- 
ware allowing for the best direct comparison of our 
HST work to brighter ground-based data. For in- 
dependent comparison, we also plot the transformed 
to the 6j-band counts from the SDSS Commission- 
ing Data (Yasuda etal. 2001). Both sets of ground- 
based counts match up well with the HST counts. 
The BBPS number count data are tabulated in Ta- 
ble 2. 

The galaxy counts for the three main morphological 
types in Fig. 3 now span a range of nearly 10 magni- 
tudes (18.5 ^ bj <J 28.5 mag), of which about 8 mag- 
nitudes (19 <; bj <; 27 mag) have decent statistics 
and reliable classifications. This is a major improve- 
ment in survey dynamic-range that could not have 
been achieved from the ground, since reliable classifi- 
cations are not possible for most galaxies in ground- 
based seeing for bj ^ 19 mag (see § 4.2). This also 
could not have been achieved from a few single deep 
HST fields, since these do not have sufficient statis- 
tics for bj ^ 24 mag. Thus, the current combined 
HST morphological counts have the potential to set 
significant constraints on galaxy evolution models in 
the flux regime 19 <; bj <; 27 mag. This is especially 



important for bj ;> 25.5 mag, where routine spectro- 
scopic measurements with 8-10 meter class telescopes 
becomes increasingly difficult and incomplete. 

The total counts are remarkably continuous and 
smooth from bj ~ 18 mag down to the formal 
HDF detection limit of bj ~ 29 mag. In the 
flux range where both surveys have good statistics 
(22 <^ bj <; 24 mag), the type-dependent counts for 
the BBPS and BBDS samples show good continu- 
ity, adding confidence that the I-band classifications 
of the faint BBDS galaxies are consistent with the 
brighter counts. 

For comparison with this new data, we plot in Fig. 3 
models for the galaxy count models that are based 
on the local LF as a function of galaxy type, as de- 
scribed in Driver etal. (1995a), Driver, Windhorst 
& Griffiths (1995) and Driver etal. (1998). These 
model have been updated using the currently favored 
flat cosmology of (Qm,^a) = (0.3,0.7) and more re- 
cent estimates of the local LF. The solid lines use the 
type-dependent LF's of Marzke et al. (1994; hereafter 
CfA) and the dashed line uses the more recent ones 
of Marzke etal. (1998; hereafter SSRS2). These are 
simple zero evolution models, which include only k- 
corrections with no explicit evolution in luminosity, 
number density or color. The models in the upper left 
panel of Fig. 3 are the sums of those of the individ- 
ual types in the other panels. No LF-normalization 
has been applied at the bright end. In fact, these 
new data show that such normalization is not needed, 
and that a global (i.e., the same factor for all types) 
LF-normalization would cause the models to overpre- 
dict the observed counts for both the E/SO's and the 
Sabc's. Such an LF-normalization has been used in 
the past to help explain the excess of fainter galaxies 
(for a discussion, see Driver, Windhorst & Griffiths 
1995; Marzke etal. 1998), but this is incompatible 
with the new data presented here. 

Therefore, the current data argue against a blan- 
ket LF-normalization with the same amplitude for 
all types. This directly results from the fact that 
the new BBPS data yield both morphology and bj- 
band magnitudes nearer to the normalization point 
(bj = 18 — 20 mag) than previous studies did. The 
adopted cosmology, which allows for more volume at 
high redshifts than an Einstein-de Sitter (Q = 1) 
model, has little effect near the normalization point, 
but serves to bring the models slightly closer to the 
data at fainter fluxes (e.g., the two models differ by 
a factor of 2 at bj ~ 26 mag). From the high dy- 
namic range of the current survey it is now becoming 
clear that the FBG excess is almost entirely due to 
the late-type galaxies, which have the steepest slope 
at the bright end, and are the most numerous type at 
the faint end of the counts. In fact, we will argue that 
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more than 90% of the FBG excess for bj ^ 25 mag is 
attributed to these late-types. The next question to 
address is how to explain this. 

4.1.1. Evolution Models vs. Renormalization 

As Fig. 3 shows, the local LF's with zero-evolution 
do not reproduce the observed galaxy counts as a 
function of type. Though this is most pronounced at 
the faint end, it is also true at the bright end, espe- 
cially for the late-types. This difference between the 
local (i.e., where the LF's have been measured) and 
the more distant Universe needs to be further inves- 
tigated. One way to model this discrepancy would 
be through a normalization of the models, which is 
the equivalent of compensating for the possibility that 
the local present day Universe is either over- or under- 
dense in a particular type (or every type) of galaxy. 
Another way is to assume some type of evolution 
(e.g., luminosity or number-density) in the galaxy 
populations at relatively low redshifts. 

Under the assumption that the data do not justify 
an LF-normalization of the models, especially for the 
early-types, we include luminosity evolution as a first 
step to explain the galaxy excess at faint magnitudes. 
As described in Driver et al. (1995a), we assume that 
galaxy luminosities evolve as L oc (1 + z) 13 , which 
is supported through the CFRS for early-mid type 
galaxies (Lilly etal. 1995b). As local references, we 
use the LF's of the SSRS2 (Marzke et al. 1998), which 
are the best available until the morphologically type- 
dependent LF's from the 2dFGRS and SDSS become 
available. Fig. 3 shows that using the CfA LF's leads 
to the same qualitative conclusions in the discussion 
that follows. The new evolutionary models are plot- 
ted in Fig. 4. The (3 = (no evolution) models are 
shown in black, (3 = +1, +2, +3, etc (positive evolu- 
tion) are shown in green and (3 = —1,-2 (negative 
evolution) are shown in red. 

The no-evolution models ((3 ~ 0) clearly provide 
reasonable fits to the BBPS data for both the E/SO 
and Sabc for bj <J 24 mag. This implies that the ma- 
jority of these types with bj <; 23 mag were in place 
by the redshift z ~ 0.5, which is the median redshift 
at this magnitude. There is an indication that models 
with a little more evolution may better describe the 
data, but a small type-dependent LF-normalization 
would also fit the data. From the bright-end to the 
limit of bj <J 23 — 24 mag, ellipticals are well fit 
by (3 ~ 1 or a modest upward LF-normalization of 
+0.1 dex. Similarly, mid- type spirals are well fit by 
(3 ~l-2 or an upward LF-normalization of ~0. 1-0.2 
dex for bj £ 24 mag. The CFRS has shown that 
giant early-mid-type galaxies underwent luminosity 
evolution with (3 ~1 since z (Lilly etal. 1995b), 
which further supports our suggestion that a signifi- 



cant upward LF-normalization of the early-mid types 
is not warranted. 

As has been noted in the past (Driver et al. 1995a), 
the late-types are clearly in excess of the no-evolution 
models over the full flux range shown here. In fact, 
even the most rapidly evolving model ((3 = +5), 
under-predicts the faint Sd/Irr counts at bj ^ 24 
mag. The best model for the total galaxy counts is 
the sum of this rapid evolution model for the late- 
types plus the no-evolution models for the earlier 
types ((3 = 0), and is shown in green in the upper 
left panel of Fig. 4. This model follows the observed 
total counts to the flux level bj <^11 — 23 mag, fainter 
than which the late-type irregulars cause the biggest 
discrepancy. If one wanted to completely fit the total 
galaxy counts with a no-evolution ([3 = 0) model, one 
would have to normalize the late-type LF by a factor 
of w 4 — 5, and such a large LF-normalization factor 
would be difficult to explain in any case as a deficit 
of galaxies in local surveys or inhomogeneties in lo- 
cal large-scale structure. A combination of luminosity 
evolution and LF-normalization is also a possibilty to 
fit the late-type counts over the full flux range pre- 
sented here, but better statistics for bj <(, 20 mag are 
required for a stronger statement to be made in this 
regard. 

In order to illustrate that the FBG excess is dom- 
inated by late-types, we consider two extreme cases. 
First, if one assumes that the zero-evolution, zero- 
renormalization models are correct, then the excess 
of the data over the models for 25 <; bj <; 26 mag 
is composed of approximately 95% Sd/Irr's. In the 
other extreme, where the E/SO's and Sabc's are mod- 
eled by zero-evolution, zero-renormalization and the 
late-types are renormalized a factor of 10 (an upper 
limit for illustrative purposes), then the FBG excess 
in the same magnitude range is only slightly reduced 
to 93% late-types. In other words, assuming that the 
models and classifications are reasonable, both the 
faint counts and the FBG excess over the model pre- 
dictions are dominated by late-types. 

One caveat that should be mentioned is that the 
morphological LF's from the SSRS2 are separated by 
type (E/S0, Sabcd, Irr/Pec) differently than our data 
where we place the Sd's in the late-type category. 
Assuming an upper-limit of Sd galaxies comprising 
25% of all spirals, this would argue for lowering the 
models for mid-types by this much and raising the 
late- type models by a similar amount. This is too 
small an effect to be seen in the current data, but 
ultimately it is an issue that should be resolved. An- 
other possible explanation for the excess of late-types 
is that they were somehow excluded from the local 
surveys, many of which were based on photographic 
data which have inherent low surface-brightness se- 
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lection issues. In fact, if one examines the data for 
the late-type LF (i.e., Marzke etal. 1998), it is seen 
that there are far fewer of this type of galaxy, rela- 
tive to the other types, in each luminosity bin. This 
means that whether the lower SB galaxies exist and 
were undetected or don't exist at all, the late-type LF 
is necessarily the least well determined of the three. 
This potentially could be resolved with a new deter- 
mination of the local LF's from more modern data 
such as SDSS and 2dFGRS, although careful atten- 
tion must be paid to not select against these late-type 
galaxies which is clearly a difficult problem to avoid. 

It is also readily apparent from Fig. 4 that none of 
the models fit beyond bj ^ 24 mag. This probably 
implies that we are seeing the epochs where galaxy 
evolution is more dramatic and not well described by 
our model. Merging and hence number density evolu- 
tion probably play a major role in this regime and our 
models are therefore most likely too simplistic. The 
issue of merging and the effects of the A are further 
discussed in § 5.4. 

The new BBPS data shown here make it very clear 
that the issues of renormalization versus evolution 
cannot be disentangled without a statistical sample of 
galaxies with morphological types that extend to even 
brighter magnitudes. The steep slope of the late-type 
counts at the brighter end (bj ^ 20 mag) indicates 
that filling this portion of parameter space should 
provide a large step forward in modelling the galaxy 
counts. Once this is done, different evolutionary sce- 
narios can be modelled and tested. Also, questions 
of merging and /or morphological evolution should be 
further investigated, but this requires the redshifts to 
be known (c/. Le Fevre etal. 1999). Unfortunately, 
the bright-end of the counts cannot be filled in by 
HST, because the surface density of bright galaxies 
is too low to efficiently use the small HST field-of- 
view. We will show in § 6 that ground based efforts 
it bj = 19—20 mag are limited by atmospheric seeing, 
but is still plausible with the large CCD surveys that 
are now being conducted, as long as the best seeing 
(FWHM ^ l'.'O) images are used. 

4.2. The Magnitude-Effective Radius Relation 

The apparent effective radii of galaxies in the BBPS 
are shown in Fig. 5 as a function of i?-band brightness 
and ANN type. Plotted is the half-light or effective 
radius (r e ), averaged over the two pass-bands, versus 
i?-magnitude. This is the radius of the galaxy con- 
taining half the light, with no assumption about the 
shape or type of profile. It is determined from in- 
tegration of the light-profile. The dark short-dashed 
lines represent the approximate detection limits for 
the exposure depths in the two surveys. The detec- 
tion limit for smaller, unresolved objects was modeled 



as a simple Gaussian source convolved with the SEx- 
tr actor convolution mask. The limit for the more ex- 
tended objects was modeled by an exponential disk. 
The fact that the object detection is based on surface- 
brightness (i.e., a minimum number of pixels above 
a specified threshold) makes it readily apparent that 
bright extended objects could be missed by automatic 
detection procedures. This is further discussed below. 
The colored lines that are almost vertical at the faint 
end give the expectations for a local galaxy of a given 
absolute magnitude, morphological type and intrinsic 
effective radius (as determined by a type-dependent 
relation between absolute magnitude and size from 
the RC3) that is redshifted back to higher z. These 
are the same models used in Odewahn etal. (1996), 
except that here we used the £Im = 0.3, f^A = 0.7 and 
H Q = 65 km/s/Mpc cosmology. 

The general trend in Fig. 5 is that the brighter 
objects appear larger, even given the plotted com- 
pleteness limits. More interestingly is that there are 
separate, although overlapping, size distributions for 
the different morphological types. This is shown in 
more detail in the histograms of Fig. 6. It is im- 
portant to emphasize that size information was not 
directly used in classifying the galaxies. As seen in 
these histograms, the general trend is that, at a given 
flux level, E/SO galaxies are the smallest in size while 
the Sd/Irr class contains the larger galaxies. The 
arrows represent the observed medians of the distri- 
butions for the BBPS (thin arrow, solid line) and for 
the BBDS (thick arrows, dashed line). In rows where 
both data sets are shown, the BBDS histograms have 
been multiplied by the ratio of the areas of the two 
data sets, while the actual BBDS numbers are shown 
on the right axes. In all cases, the BBDS medians 
are (somewhat) larger than the BBPS medians. The 
fact that the BBDS and BBPS medians do not match 
in the second and third rows is due to the fact that 
the detection limit results from the different surface- 
brightness completeness thresholds in the two types 
of surveys, which reach rather different depths. The 
BBPS limit shown in Fig. 5 shows that a galaxy with 
bj f» 24.5 mag and r e ^ Of! 8 would escape detection 
in the BBPS, but would easily be seen in the BBDS 
catalogs. Therefore, in rows 2 and 3 of Fig. 6, the 
BBPS median effective radii are lower limits, while 
the BBDS medians are more realistic, but suffer from 
small number statistics. The bottom two rows of 
the histograms come from the BBDS alone and suf- 
fer from similar SB-selection problems as the BBPS 
has at brighter levels. In general, the main result is 
that at a given brightness level, the early-type galax- 
ies appear smaller than the mid-type galaxies which 
appear smaller than the irregulars, and that fainter 
galaxies in general appear smaller than brighter ones. 
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All this can be qualitatively seen in Fig. 5 — bearing 
in mind the respective SB-completeness limits of the 
different samples — and quantitatively in Fig. 6. 

The trend for the two samples is the same to the 
completeness limit of the BBPS sample (bj <; 23.5 
mag). The two BBDS fields are also included because 
they go deeper (formally bj ^ 27.5 mag for 53W002 
and bj ^ 29.0 mag for HDF-N, see Odewahn etal. 
1996), and because there are enough galaxies past 
bj ^ 24.5 mag to make the comparison. The reason 
that the BBDS samples in Figs. 5-7 do not appear to 
go as deep in as they did in Odewahn et al. (1996) is 
that these deep survey fields were re-analyzed in an 
identical way to the BBPS, in order to make the opti- 
mal comparisons in the flux range where both surveys 
overlap. This includes the requirement that, because 
of the limited number of parallel orbits available, an 
object had to be detected in both the B and the I- 
band filters in order to be declared real. The mid- 
dle panels (23 <; bj <; 25.2 mag) for the three types 
clearly show that the larger BBPS objects are missing 
at this flux level, which is due to the stated surface 
brightness detection limits. The fact that the medi- 
ans in the third row for the BBDS are larger than for 
the BBPS in the second row clearly illustrates that 
the completeness limit is a function of size as well as 
brightness, i.e. a surface brightness limit. 

This surface brightness limit is very important be- 
cause galaxy morphology is correlated with surface 
brightness. The later-types are more extended than 
the early-types at a given magnitude, and therefore 
are of lower surface brightness. This would imply 
that the galaxy counts from the BBPS could be un- 
derestimated for the later-types for bj ^ 23 mag. 
It is not clear whether this is actually the case in 
the galaxy counts of Fig. 3. Although the BBPS 
counts are slightly lower at this flux level than in 
the deeper fields, which also have a better SB thresh- 
old, the difference is not much larger than what we 
have seen from field-to-field variations. A careful ex- 
amination of Fig. 5 shows that a small, but not in- 
significant number of bj <J 23.5 mag lower surface- 
brightness (large r e ) galaxies from the deeper fields 
would be missed in the BBPS due to the different SB- 
thresholds and these are mostly later-typed galaxies. 
Though only about a 10% effect at bj 23.5 mag, 
it quickly becomes a larger problem as one goes to 
fainter fluxes. Therefore, the surface brightness se- 
lection criteria, which was implemented by specifying 
a certain area above a threshold in SExtractor, could 
cause later-types to be omitted, but that appears to 
be a minor effect for bj <J 23.5 mag in the BBPS. 
However, this is an important issue and the surface 
brightness-morphology correlation must be consid- 
ered when detecting objects in this manner. Though 



presented here in the context of our bj ;> 19 mag 
HST data, this issue is equally important for both 
past and new local galaxy surveys, which need to be 
careful not to select against late-type galaxies. 

4.3. The Color-Magnitude Diagram 

Fig. 7 shows the (B-I) color-magnitude diagram for 
this data. The two BBDS fields were (re-)analyzed in 
exactly the same way as the 29 BBPS fields to allow 
optimal comparison. It is important to re-emphasize 
that the color was not an input parameter for the 
ANN classifier. Colors were measured in matched and 
registered elliptical apertures whose sizes and shapes 
were determined from the S-band images. For any 
given type, there is no tight correlation of color and 
magnitude that might be expected if one were ob- 
serving a uniform population at different redshifts, 
suggesting that we are seeing a scatter which is cos- 
mic in origin. 

In general, the B vs. (B — I) color-magnitude dia- 
gram shows that the reddest galaxies are E/SO's while 
the bluest ones are Sd/Irr's, with the Sabc's filling in 
the broad part middle of this range. However, there 
is also a fraction of red galaxies that are classified 
as "late-types." These could either be misclassifica- 
tions, or possibly galaxies with some star-formation 
that is reddened by dust. 

There are also some very blue galaxies that have 
been classified as ellipticals, as also noted by Driver, 
Windhorst & Griffiths (1995), who used (V - I) col- 
ors. This could result from several factors. First, 
they could be misclassified objects. This is unlikely 
at brighter flux levels, unless the profile is contami- 
nated by a neighboring galaxy, because ellipticals are 
so compact, which causes the S/N in the light profile 
to be rather good. Secondly, some of these, especially 
the brighter ones, may be misclassified and/or pos- 
sibly saturated stars. Visual inspection shows this 
to be unmistakably true only for the few bluest out- 
liers. The third possibility is that some of the "el- 
liptical" galaxies have an additional and significantly 
younger stellar population. This could be indicative 
of a burst of recent star-formation, as is seen in Com- 
pact Narrow Emission Line Galaxies (CNELGs, cf. 
Koo etal. 1995). A fourth possibility, is that these 
objects, especially those with a small measured r e , 
could have a significant AGN component. Lastly, this 
could be similar to the recent evidence against the 
traditional single burst model for star-formation in el- 
liptical galaxies (cf. Glazebrook et al. 1998, Jimenez 
etal. 1999). This would require significant ongoing 
star-formation (possibly induced by a recent minor 
merger), while the light-profile of the product had 
already sufficiently settled into an r 1 / 4 -like profile. 

To test these conjectures we have adapted the mod- 
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els of Windhorst etal. (1994), which were based on 
those of Windhorst etal. (1985) and and Kron, Koo 
& Windhorst (1985), for our filter set and cosmol- 
ogy. In each panel of Fig. 7, we have plotted these 
models for Mb = —20.7 mag which is approximately 
L* for this cosmology. The solid lines are for ages 
(i.e., last major epoch of star formation) of 14 Gyr, 
the dashed lines are for 13 Gyr and the dot-dashed 
line is a non-evolving model. The red lines are 1-Gyr 
early-burst or C-models that would represent early 
forming E/SO's that form all of their stars during the 
first Gyr after formation begins. The green lines are 
passively evolving fi = 0.7 models, where 70% of the 
mass forms into stars in the first Gyr with an expo- 
nentially declining SFR and is representative for the 
mid-type Sabc's. The blue lines are actively evolving 
fi = 0.4 models, where only 40% of the mass forms 
into stars in the first Gyr, a reasonable approximation 
for Sd/Irr's. We plot all models for each morpholog- 
ical type to allow a comparison of the morphologies 
with star formation histories. 

Note that some galaxies appear "redder" than the 
red /green upper envelope model. Another way to 
look at this is that they are instead more luminous 
than the Mb = —20.7 mag models plotted. Had one 
plotted instead the brighter Mb = —21.7 mag mod- 
els, or an even brighter model like Mb = —22.2 mag, 
these objects would also have been represented by 
these upper-envelope models (i.e., the models would 
simply shift upward/brighter by —1.0 or —1.5 mag). 
The same is true for the blue model representing the 
late-type objects — they do in general explain the ob- 
served blue objects rather well except for the few very 
blue objects seen at (B — I) <; 0.0 mag. Those would 
generally require later star formation than is present 
in the \i = 0.4 model, have a weak AGN component 
and/or possibly have larger photometric errors. 

Therefore the observed spread in color present in 
the data is roughly represented by a reasonable range 
of plotted models (the C-models though the /i = 0.4 
models), except for the very reddest and bluest ob- 
jects. The former can be explained by assuming 
M B = -21.7 mag to M B = -22.2 mag, the latter 
by assuming a more constant star formation scenario 
(perhaps fi = 0.1 — 0.2). 

5. DISCUSSION 

Here we discuss the properties of the three morpho- 
logical types as seen in the present BBPS data and 
compare it to other studies. The discussion will focus 
on the brighter BBPS galaxies, since that is the main 
subject of the present paper. 

5.1. E/SO's 

Adding the brighter elliptical counts from the 
BBPS does not substantially change the conclusion of 



Odewahn et al. (1996) that these counts can be mod- 
eled rather well by the no-evolution predictions for 
bj <^ 25 mag. The new data strengthens this conclu- 
sion and also argues against re-normalizing the mod- 
els of the E/S0 galaxy counts, although there is the 
hint that these galaxies may be undergoing some evo- 
lution at the faint end (bj ^ 24 mag). As discussed in 
§ 1, the issue of evolution of the ellipticals, which in 
principle should be the simplest to study, is far from 
settled. It is clear that the high resolution of HST is 
needed to select a sample of field ellipticals, and that 
redshifts for all objects are also needed. This was 
done for 46 galaxies from the CFRS/LDSS (Schade 
etal. 1999), which showed that field ellipticals are 
not composed entirely of old stellar populations, but 
were largely in place since z ~ 1 and therefore sup- 
port the view of an early formation epoch with occa- 
sional, more recent episodes of star-formation. 

Though the reddest galaxies at a given magnitude 
are generally of early-type, Fig. 7 shows that there is 
a broad range of colors for galaxies of elliptical mor- 
phologies. The same conclusion was reached by the 
CFRS/LDSS, who computed rest-frame (U — B) from 
the observed (V — I) using the measured redshifts. 
As previously mentioned, the bluer ellipticals could 
be similar to the CNELGs (i.e. compact and blue), 
which are thought to be distant analogs of local HII 
galaxies as well as progenitors of today's spheroidal 
galaxies (Koo etal. 1995). 

The size distribution in Fig. 6 shows the well known 
fact that the ellipticals are the most compact galaxies 
of all Hubble types. The average effective radius for 
bj <^ 24 mag is r e ~ / .'35. It was also shown in Fig. 6 
that their effective radii increase with brightness. 

Overall, these results are somewhat contrary to the 
view that all ellipticals are just old fossil galaxies, 
which all formed at the same time, whose stellar pop- 
ulations are simply passively evolving or are well un- 
derstood. While this seems true in the local Universe, 
the evidence is mounting against old red ellipticals at 
higher redshift, where there seems to be a larger color 
scatter for all types than is seen locally at the same 
rest-frame wavelengths. Even in high redshift cluster 
studies (e.g. at z = 0.83; van Dokkum etal. 2000), 
there are cluster galaxies with redder colors that are 
not of elliptical morphology. Therefore, it is danger- 
ously uncertain to use a single color for high redshift 
galaxies as an indicator of morphological type. How- 
ever, it may still be possible to determine types by 
using multiple colors which means that the galaxy 
spectrum is being sampled in more than two places 
across the SED, assuming that the spectra of the ob- 
jects are well understood. Given that there seems to 
be some discrepancy between the properties of the 
local and distant elliptical galaxies, a comprehensive 
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comparison of their SED's seems warranted in order 
to further investigate the situation. 

5.2. Sabc's 

The counts of the BBPS spirals also do not seem to 
require a large amount of luminosity evolution nor a 
renormalization to fit the models based on local LF's. 
There are far more spirals than ellipticals, but the 
numbers of mid- and late-types are about the same 
for a given magnitude. The sizes of the spirals are 
larger than the early-types and they have a median 
of r e ~ 0?41 at bj <; 24 mag. It is important to re- 
iterate that a comparison of the deeper BBDS sam- 
ple and the BBPS in a flux range where they overlap 
shows how the surface-brightness selection biases the 
sample against the larger galaxies at a given magni- 
tude. This is a result of defining the object detection 
limit as a minimum area above some threshold SB. 
This is seen in Fig. 6, where the higher signal-to-noise 
BBDS fields have a larger fraction of larger mid- and 
late-type galaxies than the BBPS for bj ^ 23 mag. 
Lilly et al. (1998) showed that the distribution of the 
physical sizes of larger (a -1 > 3.2 kpc) disk galaxies 
is roughly constant for z 1. However, a study of the 
HDF-N disk galaxy sizes showed an excess (in rela- 
tion to a CDM model prediction) of faint (Mb > — 19 
mag) disk galaxies with smaller sizes (Rj < 2 kpc) 
seen from z <J 3 even down to z ~ 0.5 (Giallongo 
et al. 2000). It is clear that the issue of size evo- 
lution and the growth of disks in spirals is far from 
settled, and the significant number of BBPS spirals 
found in the current work will allow further spectro- 
scopic studies to elucidate the physical properties of 
these objects. 

The observed colors of the spirals show a similar 
distribution to the irregulars. The spiral distribution 
is also broad, but the median is bluer than that of 
the early-types — (B — I) m ed = 2.8 mag for E/SO 
versus (B — I) m ed = 2.1 mag for Sabc for bj <; 24 
mag. This large dispersion in color is likely due to 
an intrinsic difference in the galaxies, and not en- 
tirely to bandpass-shifting. This is apparent in the 
redshift dependence of the (B — I) color given by 
Roche etal. (1997), which shows that the dispersion 
in color expected for a given type is small for z <; 1, 
especially for the later- types. A large dispersion in 
the rest-frame (U — V) colors was also shown to exist 
in the disk dominated galaxies of the CFRS/LDSS 
with only a weak dependence on redshift (Lilly et al. 
1998). 

5.3. Sd/Irr's 

The counts of galaxies beyond bj ^ 24.5 mag are 
dominated by the later- types. An examination of 
the counts of the brighter BBPS galaxies shows that 



somewhat brighter than bj ~ 22 mag, the counts be- 
come dominated by spirals. It will be very important 
to carry out this study to even brighter magnitudes 
(bj 19 mag — with HST or from the ground in 
the best possible seeing) to see exactly where this 
crossover occurs. This will be very near the region 
where the count models are normalized (bj « 18 
mag), which makes it all the more important to un- 
derstand the LF normalization as a function of galaxy 
type. 

The Sd/Irr's are generally more extended than the 
spirals and ellipticals and, by definition, have a less 
regular appearance. At bj <J 24 mag, their median 
size is r e ~ / .'61 and the SB selection effect men- 
tioned above is therefore even more important here. 
The color distribution is very similar to the spirals. 
In fact, the bright spirals seem very similar to the 
bright late-types in terms of number counts, colors 
and sizes. However, given that they are classified dif- 
ferently by the ANN, they have differently shaped az- 
imuthallty averaged light profiles. These Sd/Irr light 
profiles are similar to those of galaxies that have been 
visually classified as late-types (i.e. the training set of 
Odewahn etal. 1996). The ANN types used here are 
based on a system that does not have a merger class. 
Mergers are most effectively studied by combining the 
high resolution imaging with measured redshifts (Le 
Fevre etal. 1999) to rule out chance superpositions 
of pairs. Hence, without redshift information for all 
1800 BBPS galaxies to bj = 24.5 mag, we will de- 
fer studies of the pair-fraction and merger rate until 
spectroscopic or good photometric redshifts become 
available. 

5.4. A Possible Explanation for the Excess of All 
Types at bj ^ 24 mag 

Perhaps our most curious result is that none of our 
new models fit the observed galaxy counts func- 
tion of type at flux levels fainter than bj ;> 24 mag 
(Fig. 4). This is despite the fact that these models 
now contain the latest A-dominated cosmology and 
the best available local LF as a function of type. For 
all galaxies, there appears to be a significant excess 
for 24 <; bj <J 27 mag, which corresponds to a range 
in median redshift of approximately z me d — 0.5 to 
Zmed — 1 — 2. This excess amounts to up to a fac- 
tor 3-4 for the early-mid types, and up to a factor 
of 6-10 for the late types with respect to the local 
LF plus the best fit luminosity evolution model for 
bj <J 24 mag (i.e., (3=1 for early-mid types and 
= 5 for late types). This appears to be a robust 
result, since the model galaxy counts as a function 
of type approximately fit the brightest available data 
points, which are now around bj ~ 19 — 20 mag. 

In this section we suggest a possible explanation for 
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this, namely that we are witnessing a global and type- 
dependent excess of galaxies of all types at bj ^ 24 
mag, where the median redshift is z me d ^0.5 — with 
a larger excess for the late types — and explore pos- 
sible physical causes for this excess. 

Together with the small galaxy sizes (Figs. 5-6) 
seen at faint magnitudes, this excess suggests that 
faint galaxies may be more numerous and smaller at 
z ~ 1 — 2 compared to the ones seen today. While 
the luminous disk galaxies seen at z ^ 1 have approx- 
imately the same size distribution as that seen today 
(Lilly etal. 1998), this may no longer be true for 
z ^ 1. Hierarchical formation scenarios predict larger 
numbers of smaller objects in the epoch z — 1 — 2 and 
beyond. Based on the redshift distribution as a func- 
tion of morphological type, Driver et al. (1998) sug- 
gest that an excess of mid-late type galaxies is seen 
for I ;> 23 mag (bj ;> 24.5 mag), especially in the red- 
shift range z ~ 1 — 2. They tentatively identified this 
epoch as the epoch of disk formation. If true, then 
both the intrinsically smaller disk galaxies at z ;> 1, 
and the temporary end-products of the mergers (i.e., 
the early- types), as well as the building blocks and 
mergers-in-progress with temporarily enhanced star- 
formation (i.e., the late types/irregulars/peculiars) 
may all be enhanced in numbers past 2^1, and also 
be smaller in size and mass. 

Galaxy formation may have proceeded such that, in 
the redshift range z — 1 — 2 and at higher redshifts, 
both the early- and mid- and late-types were present, 
but likely in proportions that slowly changed with 
cosmic time, and likely with considerable migration 
back and forth between galaxy classes. For example, 
mergers between two spiral disks, or between a spiral 
and an irregular would eventually result in a bulge 
dominated galaxy, but since star-formation is known 
to not be a very efficient process, not all gas would 
be used up during that merger, nor would that gas 
necessarily reach escape velocity, so the gas would 
eventually settle back bulge dominated galaxy 
with a newly re-formed disk (e.g., Hibbard & van 
Gorkom 1996). Hence, the end-product of a merger 
would temporarily be a bulge-dominated galaxy, but 
it could grow a disk back in the next few gigayears 
after the merger, and then look like a spiral galaxy 
until the next (major) merger occurred. Similarly, 
many luminous irregulars and peculiars may be the 
temporary stages when observing mergers in action 
before a system settles as a bulge-dominated or a 
disk-dominated galaxy (Barnes & Hernquist 1996), 
although a good fraction of the late types that we 
see may likely just be the numerous smaller galaxy 
building blocks from which the hierarchical merging 
started (Pascarelle etal. 1996). In conclusion, in the 
hierarchical scenario, there would be considerable mi- 



gration back and forth between galaxy classes, and 
galaxies at high redshifts would be smaller and more 
numerous than those seen locally (Fig. 4), consistent 
with the excess seen for all types at bj ^ 24 mag 
(Fig. 4b- 4d). 

The one remaining issue that begs an explanation 
is what physical mechanism could explain the larger 
numbers of all types at bj ^ 24 mag, or z ^ 0.5-1.0? 
In hierarchical formation scenarios (Navarro, Frenk 
& White 1996, etc), bulges form relatively quickly in 
the epoch z — 3 — 5 and mostly via major mergers. 
Disks form later in the epoch z ~ 1 — 2, but more 
through the gradual (hierarchical) infall gas or mi- 
nor mergers. These predicted scenarios can be seen 
in the morphological redshift distributions of Driver 
et al. (1998). The merger rate was higher in the past 
by (1 + z) m where m ~ 2-3 (cf. Burkey etal. 1994; 
Neuschaefer etal. 1997; Le Fevre etal. 1999), but 
mostly so for z ;> 0.5-1.0. For the currently accepted 
values of Cl m ~ 0.3 and A ~ 0.7, the A-driven acceler- 
ation starts dominating the expansion of the Universe 
for the first time at z £ 1. We hypothesize that - 
as a consequence — the galaxy merger rate gradu- 
ally winds down in the epoch z ~ 0.5 — 1.0. For 
instance, groups of smaller galaxies or sub-galactic 
units that were nearly virialized at z ^ 1 will still 
virialize for z ^ 1, but groupings of such objects that 
were not even close to turn-around at z ~ 1 will be 
still expanding with the Hubble flow at z ~ 0.5, and 
probably forever do so in a A-dominated universe. 
At z <; 0.5, these late-types would never do much 
further merging, but just fade away (cf. Ferguson 
h Babul 1998). The end-result is that one observes 
a slowly evolving Universe consisting of E/SO's and 
Sabc's for z <; 0.5 (or bj <; 24 mag) — as we ob- 
serve here in Fig. 4b - 4c — plus the relatively rapid 
dwindling-away of late-types for z <^ 0.4-0.5, explain- 
ing their steep counts for b j <; 24 mag (Fig. 4d) . And 
one would observe a vast increase in numbers for all 
types at z ^ 0.5 (or bj ^ 24 mag), and especially at 
z ^ 1 (or bj ^ 26 mag), because most groupings of 
smaller objects had still plenty time to turn- around 
from the Hubble flow and overcome the effects of A 
at z ^ 1. Since the merger rate was likely much 
higher at z ;> 0.5-1.0 than at z ~ 0-0.4 (Le Fevre 
etal. 1999), merging proceeded rapidly and success- 
fully for z ;> 0.5, and vastly reduced the galaxy num- 
bers with time at z — 0.5-1.0, so that larger numbers 
of all types are seen at z ^ 0.5 — 1.0, with the largest 
increase for the late types. 

In conclusion, the new BBPS data shown here 
makes it clear that the issues of LF-normalization ver- 
sus evolution cannot be disentangled without a statis- 
tical sample of galaxies with morphological types that 
extend to even brighter magnitudes. The steep slope 
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of the late-type counts at the brighter end (bj <^ 20 
mag) indicates that filling this portion of parameter 
space should provide a large step forward in modeling 
the galaxy counts. Once this is done, different evo- 
lutionary scenarios can be modeled and tested. Also, 
questions of merging and/or morphological evolution 
should be further investigated with better statistics 
from larger surveys, with wider dynamic range, and 
through a systematic assessment of the effects from 
the uncertain rest-frame UV (Windhorst et al. 2002) 
on the classifications at the faintest magnitudes. 

6. SUMMARY AND FUTURE WORK 

We have presented an HST survey that connects 
the extremely deep HST studies, such as the HDF-N 
(24 £ bj <^ 29 mag) and existing ground based stud- 
ies such as the RC3 (bj ^ 17 mag). For example, the 
Stromlo-APM Redshift Survey (Loveday et al. 1996) 
provided a catalog with morphological types for 1797 
galaxies and is complete to bj = 17.15 mag, although 
there is some question as to the reliability of these cat- 
alogs (Pozzetti, Bruzual & Zamorani 1996). Unfortu- 
nately, there is a relatively small amount of good mor- 
phological data available in the literature for a mag- 
nitude range of 17.15 < bj < 19 mag. The bright-end 
of the galaxy counts cannot be filled in by HST be- 
cause the surface density of bright galaxies is too low 
to efficiently use the small HST field-of-view. Our 
expectation is that this flux range can be addressed 
from the ground with existing telescopes and larger 
area detectors in good seeing. The SDSS and other 
wide angle ground-based CCD surveys should soon 
provide this data. The key problem with these sur- 
veys will be in classifying the fainter (bj ;> 19 mag) 
galaxies from the ground due to seeing-related effects, 
since their median r e values are <; 1" (see Fig. 5 here) , 
and rapidly decrease towards fainter fluxes. Hence, 
the flux range 19 <^ bj <; 24 mag must be studied with 
HST, which was the purpose of the current study. 

Another important piece of missing information is 
the redshifts of the BBPS. Since the BBPS galax- 
ies are brighter than bj <j 25 mag, this is a project 
that can be started on a 4 meter class telescope and 
finished on an 8-10 meter class telescope. Objects 
of known redshift (even if estimated through photo- 
metric redshifts) can, in principle, be more accurately 
classified, because one can more effectively correct for 
the effects of band-pass shifting, mentioned in § 3.5. 
The BBPS along with measured redshifts can pro- 
vide a wealth of information when combined with 
the wealth of HDF-N redshifts which are now mea- 
sured (Cohen etal. 2000). Given the photometric 
redshift distributions of Driver etal. (1998) for the 
HDF-N, it is clear that the majority of the galaxies 
at bj <; 25 mag studied in this paper are at z <; 1. 



Therefore, in order to classify galaxies observed in the 
5-band, we need to know what local galaxies look 
like in the rest-frame [/-band or at slightly shorter 
wavelengths (2500-3000^4). Relatively recent stud- 
ies of the near-UV morphology of nearby galaxies 
show some morphological differences as compared to 
the £>-band (Giavalisco et al. 1996; Burg et al. 1997; 
Kuchinski etal. 2001; Marcum etal. 2001; Wind- 
horst etal. 2002), but a full quantitative analysis of 
how this would affect the galaxy counts has yet to 
be performed. The /-band classifications used in the 
present paper largely avoid this issue, except for the 
relatively small number of higher redshift {z ^ 1.5) 
galaxies whose classifications are therefore necessarily 
uncertain. 

The use of artificial neural networks was previously 
shown to be effective in classifying a large number of 
galaxies in a quantitative and systematic (i.e., repro- 
ducible) way (Odewahn etal. 1996). This method, 
based almost solely on the shape of the measured 
light-profiles, has been applied to our much larger 
data set at brighter levels. The appeal of this method 
is that it uses a large number of photometric param- 
eters, and that it is also "trained" based on human 
classifiers in an effort to categorize the actual appear- 
ance in a systematic, albeit non-linear way. Other au- 
tomated techniques in use today can produce consis- 
tent results and provide other useful types of informa- 
tion, while also providing a good consistency check on 
the method used here. The next logical step is to see 
if the ANN method can be improved by using some 
of the 2-dimensional information in the images. The 
ultimate goal, which no published method to date has 
achieved for faint galaxies (including the one used in 
this paper), would be to quantitatively measure the 
true morphology in an automated way. This would 
involve distinguishing between more subtle features 
such as spiral arms, bars and rings, as well as differ- 
entiating between, for example, Sa and Sb, and trac- 
ing their behavior with redshift (see Odewahn et al. 
2002). 

In summary, the galaxy counts, size distribution, 
and (B — I) color distributions seen in the deep HST 
studies are consistent with what we are now seeing 
with good statistics for the brighter BBPS galaxies. 
The excess of faint galaxies for bj ^22-23 mag is 
dominated by the late-types (cf. Driver et al. 1995a, 
Driver, Windhorst &: Griffiths 1995). There are rel- 
atively few early-type galaxies at faint magnitudes. 
Models indicate that either luminosity evolution or 
an extra dwarf population of late-types is needed to 
explain the counts of the later types. Redshift surveys 
suggest the former, i.e. evolution through episodic 
starbursts. The new data and new models presented 
here do not support the need for re-normalizing the 
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total galaxy count models at bj ~ 18 mag. While 
not ruling out this need, they show that if this re- 
normalization is necessary then it must be a function 
of galaxy morphology. Brighter objects appear larger 
with the broad trend of increasing apparent size as 
one goes from early to mid- to late- types. This ap- 
pears to be true over a range of almost 10 magnitudes. 
There is no sharp size cutoff between types. In gen- 
eral, the early-types are redder, while the later types 
are bluer with the mid-types in between. Again, there 
is no simple way to differentiate types based on ob- 
served colors, especially with only two bands. 

From the current study, we provided the first sys- 
tematic 6j-band counts as a function of galaxy type 
to address the problem of normalizing the model 
galaxy counts, which use the known local LF's as a 
function of morphological type. The galaxy statis- 
tics at the bright end (bj <^ 19 mag) are still rather 
poor. At bj ~ 18 mag, the galaxy counts are approxi- 
mately 100 galaxies per 0.5 mag per square degree. A 
single WFPC2 field is 0.0013 square degrees, which 
means that about ten fields are needed to see even 
one galaxy at random in this magnitude bin. This 
implies that a few hundred fields would be needed to 
have adequate bright end statistics. Therefore, larger 
detectors are needed. We expect that we can 
study the magnitude range 16 bj <^ 19 mag from 
the ground (S. H. Cohen, in preparation), but only 
in good seeing, as e.g., using images from the NOAO 
Deep Wide-Field Survey (Jannuzi & Dey 1999), and 
soon also from the Millenium Galaxy Catalog (Liske 
et al. 2002), when classifications are added to it. This 
may be feasible because the effective radii at this 
magnitude are larger than the seeing disk in good see- 
ing conditions (see § 4.2). The exact magnitude limit 
to which these classifications can be reliably pushed 
from the ground is not yet known, but our expecta- 
tion is to get complete and reliable classifications to 
bj ^ 18 — 19 mag. This combined data set will pro- 
vide better statistics at the bright end, which then 
can be used to more firmly address the normalization 
problem at bj = 18.0 mag discussed earlier, and de- 
termine if the normalization factor is a function of 
galaxy type. 

Some of the requirements to improve upon the in- 
terpetation of the faint galaxy counts in order to truly 
get a handle on issues of galaxy formation and evo- 
lution are as follows: 

1. Brighter galaxy counts (16 <; bj <^ 19 — 20 mag) 
as a function of morphological type 

2. A better handle on the morphological classifi- 
cation accuracy 

3. Better statistics throughout 



4. CCD-based type-dependent LF's classified us- 
ing consistent methods 

5. Consideration of surface brightness selection ef- 
fects 

6. Measured redshifts and better determined k- 
corrections 

An important point to close on is that we have an 
incomplete understanding of the local and interme- 
diate distance Universe. Our knowledge and inter- 
pretation of the distant high-z Universe will always 
be limited by this. The advent of bigger and bet- 
ter telescopes brings about the temptation to observe 
the fainter and more distant objects in the Universe. 
These studies still need to be complemented by those 
of more nearby objects, such as that which is pre- 
sented here, in order to paint the complete picture. 
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Fig. 1. — /-band magnitude vs. effective radius for the BBPS. This figure shows that misclassification between stars and 
galaxies occurs for I > 22 mag. At this magnitude, a fraction of stars may be classified as early-type galaxies. Note that stars are 
only plotted for the BBPS fields and not for the two BBDS fields. The dashed lines indicate the approximate completeness limits 
for the BBPS and BBDS, as discussed in the text. The objects with pluses superposed were re-classified as "stars", and are likely 
globular clusters in the Virgo cluster, as explained in § 3.6. The scale indicated across the top is for r e in arcseconds. 
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Fig. 2. — Differential star counts for the BBPS data. Objects initially classified as stars are shown as upward pointing arrows 
and the final star counts are shown as stars with error-bars (see § 3.6). The solid line shows the predicted star counts from the 
Bahcall & Soneira (1981) Galaxy model averaged over the Galactic coordinates of all 29 different HST parallel pointings. For 
comparison, the total counts for objects classified as galaxies (triangles) and those classified as E/SO's (circles) are also shown (see 
also Fig. 3). Note that to bj < 23.0 mag, the initial star-galaxy separation is good. Beyond this limit, the galaxies begin to far 
outnumber the stars, so that a few misclassifications of stars will not substantially affect the results on faint galaxies. However, 
at bj > 23 mag, the star counts become less reliable. 
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Fig. 3. — Differential B-band galaxy counts for all BBPS galaxies as a function of morphological type and bj magnitude. The 
solid blue points represent the new BBPS data that cover 18 <C bj <J 24 mag. The open triangles (HDF-N) and open circles 
(53W002) represent the counts from the two deeper fields analyzed in the same way as the BBPS counts. The curves are predictions 
from local LF plus no-evolution models as described in Odewahn etal. (1996) and in Driver etal. (1995a). The models assume 
a cosmology with (Qm,Qa) equal to (0.3,0.7). These non-evolving models use the type-dependent LF's of Marzke etal. (1994; 
dashed lines) or Marzke etal. (1998; solid lines). The open squares show the upper limits at the bright-end, as described in the 
text. The red diamonds are the total ground-based CCD counts from the Millenium Galaxy Catalog (Liske et al. 2002) and the 
green squares are from a first part of the SDSS (Yasuda etal. 2001). These SDSS and MGC total counts are perfectly consistent. 
The reliable classification limit in the BBDS is bj <C 27 mag, based upon the available training-sets and S/N in the images (see 
text and Odewahn etal. 1996). 
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Fig. 4. — Evolutionary models for the differential S-band galaxy counts as a function of morphological type and bj magnitude. 
The solid blue points represent the new BBPS data that cover 18 < bj < 24 mag. All plotted data symbols are the same as 

in Fig. 3. These models use the LF's of Marzke etal. (1998), and assume galaxy evolution in the form of L oc (1 + z) 13 . The 
no-evolution (f3 = 0) models are shown as solid black lines, positively evolving models with /3 = +1, +2, +3, +4, +5 as green lines, 
and negatively evolving model with /3 = —1, —2 as red lines. In panel (a), the red line is for /3 = for all types, and the green line 
is for /3 = for the E/SO and Sabc galaxies and j3 = +5 for the Sd/Irr population. Clearly the Sd/Irr counts in panel (d) are well 
above any reasonable model for bj > 23 mag, and this is causing most of the excess of the total counts in panel (a). 
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Fig. 5. — B-band magnitude-effective radius relation for galaxies in BBPS and BBDS data sets. The thick dashed lines 
represent the approximate completeness limits for the deep and shallower surveys as described in the text. The lines that are 
almost vertical at the faint end from Odewahn etal. (1996) show the expected fej-r e relation for redshifted RC3 galaxies of a 
given absolute magnitude and assumed cosmology, as indicated. Objects below the horizontal line at b.j = 27 mag are beyond the 
reliable classification limit of the BBDS. The scale across the top indicates r e in arcseconds. The quantity r e plotted here is the 
average of the individual effective radii measured from the light profiles in the B and I pass-bands. The r e values are in general 
very similar between the B and I filters, so yielding a better measure of r e . 
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Fig. 6. — Distributions of galaxy sizes as a function of B-band brightness and /-band morphological type. The arrows indicate 
the median effective radii for the given distributions. The solid histograms and thin arrows are for the BBPS (using the number 
scale plotted on the left), and the dotted lines and thick arrows are the BBDS (using the scale on the right). In the second and 
third rows (23 < bj < 25.2 mag), the BBDS data has been scaled up by the ratio of the areas of the two data-sets for comparison 
purposes. The completeness of the histograms as a function of flux and size is discussed in the text. 
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Fig. 7. — The B vs. (B — I) color magnitude diagram as a function of galaxy type. Each of the three different morphological 
ANN types are plotted in different panels as indicated. Both the two BBDS (solid symbols) and 31 BBPS (open symbols) fields 
were measured in exactly the same way for consistency. The slanted and horizontal black lines indicate the B-band and /-band 
50% detection limits, respectively. The solid black line at bj — 27 mag indicates the /-band classification limit (as limited by the 
uncertain rest-frame near-UV morphology and S/N considerations). The models are plotted for a Mb = —20.7 mag galaxy that 
is either 14 Gyr old (solid) or 13 Gyr old (dashed). The non-evolving model is plotted as a red dot-dashed line. The red, green 
and blue lines are meant to be representative of E/S0, Sabc and Sd/lrr galaxies, respectively, although they are not neccessarily 
unique. More details of the models can be found in the text. 
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Table 1 
BBPS Data Summary. 
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24.5 


24.75 


24.25 



1 These are the field names listed in the order the data was received. The original data can be obtained from the HST Archive 
(http://archive.stsci.edu/hst/search.php) by entering the coordinates in the search form. 

2 Galactic absorption As (in mag) from Burstcin & Heiles (1982). 

3 Total integration times in seconds (number of exposures). 

4 There is a total of about 151 Parallel HST orbits in this data set. 

5 SB(6j,lim) is the limiting surface brightness in bj mag arcsec~ 2 for a given field averaged over the 3 WF detectors. The PC limit is 
approximately 1.5 mag arcsec~ 2 brighter. 

6 Center of faintest complete 0.5- mag bin in the total galaxy counts (50% complete). 

7 Field is in or near the Coma or Virgo superclusters (see §3.6 and §3.7 for details). 
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Table 2 

Differential 6j-band Galaxy Counts as a Function of Type 
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24.0 


24.5 


4.006 


2.857 


3.556 
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25.0 


4.162 


3.021 
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1 All counts are log 10 (Number/sq. degree/0.5 mag) 

2 Lower and upper 6j-magnitude bounds of given bin 

3 Area is 0.0370 square degrees, except in last 2 bins, where 
the effective areas are 0.0361 sq. deg. and 0.0295 sq. deg., 
respectively, as discussed in §3.7 

4 The 23 <J bj <J 23.5 mag bin is the faintest bin where the 
total counts are better than 90% complete. 



